Evolution of a coherent array of Bose-Einstein Condensates in a magnetic trap 
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We investigate the evolution process of the interference pattern for a coherent array of Bose- 
Emstem condensates in a magnetic trap after the optical lattices are switched off. It is shown that 
there is a decay and revival of the density oscillation for the condensates confined in the magnetic 
trap. We find that, due to the confinement of the magnetic trap, the interference effect is much 
stronger than that of the experiment induced by Pedri et al. (Phys. Rev. Lett, 87, 220401), where 
the magnetic trap is switched off too. The interaction correction to the interference effect is also 
discussed for the density distribution of the central peak. 
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I. INTRODUCTION 

The development of the technologies of the laser trapping and evaporative cooling has yielded intriguing Bose- 
Einstein condensates (BECs) [jD, a state of matter in which many atoms are in the same quantum mechanical state. 
The realization of BECs has made remarkable theoretical and experimental advances on this exotic quantum system 
P,|J . Recently, optical lattices are used to investigate further the unique character of the ultra-cold atoms Q . The 
application of the optical lattices to the ultra-cold atoms is very promising, such as the quantum computing scheme 
proposed in [^]j6j. Experimentally, the macroscopic quantum interference effect R and thermodynamic properties 
H of the BECs in the optical lattices have been investigated thoroughly. In addition, the superfluid and dissipative 
dynamics p|JlO|| of a BEC in the optical lattices are investigated in the experiments too. In particular, through the 
application of the optical lattices, the quantum phase transition from a superfluid to a Mott insulator in Bose-Einstein 
condensed gases is observed in a recent experiment 

Recently, the expansion of a coherent array of BECs is carried out to illustrate the interference effect in the 
experiment by Pedri et al. fll2fl, when both the magnetic trap and optical lattices are switched off. In another 
experiment by Morsch et al. |13[ | , the expansion of the condensates is investigated when only the magnetic trap is 
switched off. In this paper, we give a theoretical investigation on the expansion of the condensates when only the 
optical lattices are switched off. Due to the confinement of the harmonic potential, we find that there are several 
interested phenomena not encountered in the experiment by Pedri et al. ]T^ | and Morsch et al. In the presence 
of the harmonic potential, the interference effect would be much stronger than the case when the magnetic trap 
is switched off too. In this situation, researches show that the interaction between atoms would give important 
correction, in contrast to the case investigated by Pedri et al. [ fl2|| . 

II. WAVE FUNCTION IN HARMONIC POTENTIAL AFTER OPTICAL LATTICES ARE SWITCHED 

OFF 

In the experiment conducted by Pedri et al. Il2| , the external potential of the Bose gas is given by |12| : 



V = \m (c 2 * 2 + «i ( y 2 + z 2 )) + sE R cos 2 + ~) , 



The last term represents the external potential due to the presence of the optical lattices. In the above expression, u) x 
and ujj_ are the axial and radial frequencies of the harmonic potential, respectively. In addition, A is the wavelength 
of the retroreflected laser beam, and sEr denotes the depth of the optical lattices. For the optical lattices created by 
the retroreflected laser beam, the last term has a period d — A/2. In other words, d can be regarded as the distance 
between two neighboring wells induced by the optical lattices. In this paper, the experimental parameters in fll2| are 
used to calculate various physical quantities, thus it would be useful to give them here. The experimental parameters 
in pi] are u x = 2tt X 9 Hz, oj±_ = 2tt x 92 Hz, A = 795 nm, E R = 2nh x 3.6 kHz, and s = 5. 
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Due to the presence of the optical lattices, there are an array of condensates formed in the combined potential, 
when the temperature is lower than the critical temperature. In this work, we shall investigate the case of the strong 
tunneling between neighboring BECs, which holds in the experiment by Pedri et al. In this situation, all the 

condensates are fully coherent, and can be described by a single order parameter. To emphasis the role of the optical 
lattices, our researches are carried out mainly on the character of the coherent array of condensates in the a;— direction. 

In the presence of a magnetic trap, the number of atoms in each well should be different. Based on the analysis of 
the 3D model of the condensates in the combined potential jip], the ratio between the number of condensed atoms 

in fc— th and central wells is given by N^/Nq = (l — k 2 jk\.^) , where 2&m + 1 represents the total number of 
condensates induced by the optical lattices. In this situation, using the Gaussian approximation in the a;— direction 
for each well, the normalized wave function in coordinate space takes the form 



kM ( k 2 
(fo(x) = A n ^2 ( 1 _ — ) ex P[~ ( x ~ kd f / 2<j2 }> 



fc=-fej 



(2) 



where a denotes the width of the condensate in each well. It can be calculated by numerical minimization of the 
energy of the condensates [|l2|. In the following calculations in this paper, a — 0.25c? (l2) is used for s = 5. In the 
above expression, the normalized constant A n takes the form 



A,, 
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V(16fc^ - 1) /15fc^^ 1 /4 (7 l/2 



(3) 



It is well known that once the wave function at an initial time is known, the wave function at a later time can be 
obtained through the following integral equation Jl4[ : 



cp Q (x,t) 



K(x, t;y,t = 0)<p a (y, t = 0) dy, 



(4) 



where ipo (y, t = 0) is the wave function at the initial time t = which is given by Eq. (|^), and K (x, t; yA — 0) is the 
well known propagator. For the atoms in the harmonic potential, the propagator can be expressed as pM: 



K(x,t;y,t = 0) 



rrioj , 



2n iTi sin uo x t 



1/2 



exp 



2Ti sinuj x t 



[(x 2 + y 2 ) cos uj x t — 2xy] 



(5) 



From the formulas (Q), (Q) and (||), after a straightforward calculation, one obtains the analytical result of the wave 
function confined in the magnetic trap 



t] = sin^i(ctU + i7 ) t, ( l ~ £ I eX P 



(kd cos u> x t — x) 



2a 2 sin 2 uj x t (ctg 2 w x t + 7 2 ) 



exp 



% (kd cos uj x t — x) ctguj x t 
2ja 2 sin 2 oj x t (ctg 2 uj x t + 7 2 ) 



exp 



i (x 2 cosuj x t + k 2 d 2 cosu> x t — 2xkd) 
27 a 2 shiujxt 



(6) 



where we have introduced a dimensionless parameter 7 = h/mu> x a 2 . Assuming N denotes the total number of particles 
in the condensates, the density distribution in a;— direction is n (x, t) — N \<~po(x, t)\ . 



III. PERIODICITY OF THE DENSITY DISTRIBUTION AND THE MOTION OF N = ±1 PEAK 



Due to the confinement of the harmonic potential, the density distribution in x— direction should exhibit a periodic 
character. From Eq. (|^), it is easy to find that the period of the density distribution n (x, t) is determined by uj x T = tt. 
For the experiment in p2| , lu x = 2n x 9 Hz, this means that the period of the density distribution is given by 500/9 
ms. 

In Fig. 1, displayed is n (x = 0, t) when only the optical lattices are switched off. The periodicity of the density is 
clearly shown in the figure and in agreement with the analytical result given by T — tt/uj x . We see that the density 
at x = reaches a maximum value at time t m = (2m— 1)tt/2uj x , with m positive integer. At time t m , the wave 
function takes the form 
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<Po ( 



The maximum density at x — is then 



x,t m ) = A n J — [i-TT exp 



2a 2 7 2 



{ — 2 



(7) 



n (x = 0, t m ) 
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In the case of fc^/ >> 1, the above expression can be approximated as: 

n(x = 0,t m ) Na>l_ ideal , 

where a 2 _ ldea , is given by 



M 



ideal 



(8) 



(9) 



(10) 



In the experiments, the depth of the optical lattices can be changed through the variation of the parameter s. From 
JlH , we know <r oc l/s 1 / 4 and Um oc 1/ct 1 / 5 . Therefore, n(x = 0, i m ) cx 1/s 1 / 5 . This shows that in the non-interacting 
model, the maximum density at x — decreases with the increase of the depth of the optical lattices. 

In Figures 2 (a)-(d), we show the evolution of the density distribution of the condensates confined in the magnetic 
trap, after the optical lattices are switched off. The density distributions are shown at t = 0, 0.1tt/uj x , 0.37r/o; x , and 
0. 5tt/lo x in these figures. The motion of the n = ±1 peaks is clearly shown in these figures. The oscillation motion of 
the n = ±1 peaks is due to the confinement of the harmonic potential. In fact, the motion of the n = ±1 peaks can 
be described very well using the classical harmonic motion. Using the classical harmonic motion, the motion of the 
n = ±1 peaks is determined by the following expression: 



2irh 

x n= ±i (t) = ± cos 

mw x d 



(11) 



When obtaining the above formula, we have used the fact that the momentum distribution is characterized by sharp 
peaks at the values p x — n2irh/d JL2|. The solid line in Fig. 3 shows the harmonic motion of the n = 1 peak using 
the above formula, while the circles show the result given by Eq. (pf ). We see that the classical harmonic motion 
agrees quite well with the result given by Eq. (||). In a sense, Eq. ( jll| ) describes the motion of the center of mass of 
the n = ±1 peak. Thus, we anticipate that the interaction between atoms will not change the motion of the n = ±1, 
although it will affect the density and width of the n = ±1 peak. 

From Fig. 2(d) we see that at time t m the density distribution in x— direction exhibits a very sharp peak at the 
center of the magnetic trap. The maximum density of the central peak is given by Eq. (|9|). As a comparison, assume 
there are N atoms confined in an identical magnetic trap, but there are no optical lattices to induce the interference 
effect. In this situation, in a:— direction, the density distribution at x = is given by 



0) 



2 

gmLU 2 ^ ' 



(12) 



where fJ-mag is the chemical potential of the Bose gas confined in the magnetic trap. The ratio between n (x — 0,t m ) 
and n mag (x — 0) is then 



n{x = 0,t m ) = N gmuj±a^ ideal 



2 



(13) 



For the experimental parameters in n [x = 0, t m ) /n. mag (x = 0) = 29.6. This shows clearly that there is a very 
strong interference effect for the case considered here. 

If both the magnetic trap and optical lattices are switched off, there is also a sort of interference effect. The 
maximum value of the density distribution in this case is given by 



n bs (x = 0) 



N 



(14) 
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From Eq. (|l4|), rib s (x = 0) oc s 1 / 4 . In contrast to the case where only the optical lattices are switched off, ribs (x = 0) 
decreases with the decrease of the depth of the optical lattices. 
The ratio between n^s (x = 0) and n mag (x = 0) is then 

n bs (x = 0) _ Nogmuj^ 

n m ag {% = 0) 7T 3 / 2 (T^^ ag ' 

For the experimental parameters in p^ |, we have nb s (x — 0) /n mag (x — 0) = 2.2. From Eqs. (Q) and (fil|), 
n (x = 0, t m ) /nts (x = 0) = 13.7. Therefore, when only the optical lattices are switched off, the interference ef- 
fect would be much stronger than the case when the magnetic trap is switched off too. When only the optical lattices 
are switched off, at t m the density of the central peak is very high, and we anticipate that in this situation the 
interaction between atoms would give important correction. 



IV. DECAY AND REVIVAL OF THE DENSITY OSCILLATION 



From Fig. 1, we sec that there is a phenomenon of decay and revival of the density oscillation at x = 0. We now turn 
to discuss this unique character. To proceed, it is useful to introduce an important time scale which determines when 
the interference between two neighboring condensates begins to occur. Before the magnetic trap and optical lattices 
are switched off, from the Gaussian approximation of the condensates in each well, the width of the condensates in 
each well is given by 

^. /ffi^a* tie, 

J exp [—x z /a z \ ax 

It is easy to get Axo = from the above formula. When the optical lattices are switched off, for time much 

smaller than tt/lo x , the condensates can be approximated as a free expansion and the width of the condensates would 
increase in this situation. Based on the analysis of the spreading of the wave packet, the width of each condensate is 
given by 



Ax(t) = Ax 0] jl + ^ o (17) 

When Ax (t) = d, the condensates in neighboring wells begin to interfere with each other. By setting Ax (t) — d in the 
above formula, we obtains a time scale t w which determines when the interference between neighboring condensates 
begins to occur. From Eq. (|l7|), it is easy to find that the following analytical result of t w can give a rather well 
approximation 

a dm , . 

t, - (18) 



From Eqs. ( |l7j ) and (|18|), for t > t w , one gets the following useful result: 

Ax(t) = —d. (19) 

As illustrated in Fig. 1, the oscillation of the density at x — will cease ultimately when t > kMt w - However, 
when the time approaches n /lo x , the density oscillation will reappear. Note that the time for the revival of the density 
oscillation is determined solely by the axial frequency of the harmonic potential. This shows that the confinement of 
the harmonic potential plays a crucial role for the revival phenomenon of the density oscillation. To verify further 
the decay and revival of the density oscillation, Fig. 4 shows the density at x — Analogous decay and revival 

of the density oscillation are illustrated clearly in the figure. However, the oscillation of the density disappears at a 
longer time t = 1.47/cm^, in comparison with the case at x = 0. 

In Fig. 5, we display the time of the disappearance of the density oscillation for different locations in the region 
< t < tt/uj x . We can give a rather simple interpretation for this result. When the optical lattices are switched off, 
the width of the expanding condensates in each well will increase. For the location at x = 0, there are more and 
more expanding BECs interfere at this point with the development of the time. This is the reason why the density 
at the point x = will oscillate intensely. When t > kMt w , however, all expanding BECs have participated in the 
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interference at the point x — 0. Therefore, the oscillation of the density at x = will cease at time longer than kutw- 
Generalizing this result, the time for the disappearance of the density oscillation for different locations is given by the 
following simple expression: 

t= (k M + x/d)t w . (20) 

The solid line in Fig. 5 displays the above analytical result. We see from Fig. 5 that this simple expression agrees 
well with the result given by Eq. (|^). Maybe the slight difference from the result given by Eq. (^J) lies in the fact that 
we do not account for the effect of non-uniform atom distribution in each well when obtaining Eq. (|2p|). 

It is worth pointing out that the disappearance of the density oscillation do not mean the disappearance of the 
interference effect between the expanding condensates. When < t << (2km + l)t w , there are only several expanding 
condensates interfering with each other, and there are interference fringes (or density oscillation) in this situation. For 
(2k m + 1) t w < t < tt/lu x — (2fcjvf + 1) t w , however, all expanding condensates will interfere with each other, and this 
means the emergence of the diffraction fringes. In fact, n = and n — ±1 peaks in Fig. 2 (a)-(d) should be regarded 
as the diffraction fringes, rather than the interference fringes. Note that the phenomena of diffraction and interference 
are basically equivalent. Different from the interference phenomenon, however, the diffraction phenomenon should be 
regarded as a consequence of interference from many coherent wave sources. In a sense, Eq. (^0|) gives the time for 
the emergence of the diffraction fringes, which means the disappearance of the density oscillation. 



V. INTERACTION CORRECTION TO THE CENTRAL PEAK AT T M 

In the case of non-interacting model, we have shown that, at time t m , there would be a sharp central peak in the 
magnetic trap. In this case, the interaction between atoms can not be simply omitted, in contrast to the case when the 
magnetic trap is switched off too. At time t m , the central peak of the density in ir— direction can be approximated as 
a Gaussian distribution. After the optical lattices are switched off, using Thomas-Fermi approximation in the radial 
direction, the square of the modulus of the 3D wave function at time t m takes the form 



\ip (x,r±,t m )\* 



OCc a ±_ ex P 



2x 21 
#1 



f 



R\ 



(21) 



where R\ 



y/2no/ muJ '± with /i = muj 2 k 2 M d 2 /2 In the above expression, a x — (2 / V) 1 ' 14 / '\fRx~ and a± = 

y/2/irR 2 L are the normalized constants in the axial and radial directions, respectively. Obviously, Na 2 represents 
the density n [x = 0, t m ) in the a;— direction. Due to the repulsive interaction between atoms, we anticipate that 



2,2 

a_ < a. 



x — ideal' 

Assume E int and E^ in are the interaction energy and kinetic energy of the central peak at t m , respectively. The 
interaction energy of the central peak is given by 



E u 



gN 2 



\ipa(x,rj_,t m )\ dV 



V2~9N 2 al 
3ttR 2 , 



(22) 



Assuming the total energy of the condensates is E a u , we have 

Ekin + Ei n t + Eho = E a u , 



(23) 



where Eho is the potential energy of the condensates. For the central peak, Eho can be omitted safely. In the case 
of non-interacting model, E a u will transform fully to the kinetic energy and potential energy of the condensates once 
the optical lattices are switched off. Due to the presence of the repulsive interaction between atoms, the maximum 
density of the central peak would be smaller than the result of the non-interacting model. 

Note that the kinetic energy of the central peak can not be calculated through N j ^ \^ \J\fo (x,r±,t m )\j dV, 



because the phase factor is different for different well (See Eq. ([?])). From the uncertainty relation, assume E^ 
^/R 2 oc a 4 . In the presence of repulsive interaction, we have 



Eki 



x — ideal 



-E, 



u- 



(24) 



From Eqs. (J22j), (E3T) and (E4), one obtains the following equation to determine a 2 : 
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(3 2 + 0(3 - 1 = 0, 



(25) 



where (3 = (ct x /&x-ideai) 2 ■ The value of (3 reflects how the repulsive interaction between atoms reduces the density 
of the central peak. In the above expression, the dimensionless parameter 9 — Ei nt (a x — a x -ideai) /E a u. From Eq. 
(^B|), we have 



(26) 



Now let us turn to discuss the total energy of the condensates which is necessary to calculate the value of (3. The 
total energy of the condensates can be obtained through the sum of the energy of the condensates in each well before 
the optical lattices are switched off. Before the optical lattices are switched off, the normalized wave function in k— th 
well takes the form 



fak (x, rj_) = (p ok (x) ip k (rj_) 



where 



<fok (x) 



1 

/7R7 



1/2 



exp 



(x — kdy 



(27) 



(28) 



and ipok i r ±) is the normalized wave function in the radial direction. From Eq. (p7|), we have 



E a it = ^ < N k / ^ ok (x, n 

7 L. \ J 



k— — k\ 



2m 



1 2 2,1 22 

-mui„„x + -mu 1 r 1 

2 2 



Vofe (x,r±) dV 



9NI 



\<Pok (x,rx)\ dV 



(29) 



where uj xe is the effective harmonic frequency in the x— direction of the well induced by the optical lattices. The last 
term in the above expression represents the interaction energy of the condensates in each well, and it is easy to verify 
that it can be omitted safely. In addition, the kinetic energy and potential energy in the radial direction can be also 
omitted because u xe >> ui±, and huj xe >> fiQ. In this case, one gets 



E a l, 



N 



\Ama 2 



1 22 



(30) 



For the experimental parameters in |l2j, from the formulas (22), (p6|), and (pO|), the calculation shows that (3 = 0.80. 
This shows clearly that the repulsive interaction between atoms would reduce the density of the central peak at t m . 

With the variation of the parameter s,, the maximum density of the central peak can be also calculated based on 
the method given here. In the interacting model, the solid line in Fig. 6 shows the ratio between n (x = 0, t m , s) and 
n (x = 0, t m , s = 5). The result of the non- interacting model is also shown in the figure. 



VI. DISCUSSION AND CONCLUSION 



In brief, the evolution process of the condensates is investigated after the optical lattices are switched off. We find 
that the density oscillation exhibits a phenomenon of decay and revival, based on the numerical result of the evolution 
of the density distribution. The decay of the density oscillation is interpreted as the emergence of the diffraction 
phenomenon, which is regarded as a consequence of interference from a lot of coherent expanding condensates. Due 
to the confinement of the harmonic potential, there is a periodic character of the density distribution, and it is this 
periodic character which leads to the revival of the density oscillation. In contrast to the condensates in the magnetic 
trap, there is no revival of the density oscillation, when both the magnetic and optical lattices are switched off. In 
addition, in the case of non-interacting model, it is shown that the maximum value of the density distribution at x = 
would be approximately 30 times larger than the case when there are no optical lattices to induce the interference 
effect. 

It is shown here that the repulsive interaction between atoms has an effect of reducing the maximum density of the 
central peak. In a real experiment in future, maybe the experimental result of the maximum density of the central 
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peak would be smaller than the theoretical prediction given here, because it is possible that there is a loss of the total 
energy of the condensates during the process of the removal of the optical lattices. For the attractive interaction such 
as Li, the role of interaction would become very important when the atoms are confined by the combined potential. 
For example, when only the optical lattices are removed, based on the non-interacting model, the density at x = 
would increase largely due to the interference and confinement of the magnetic trap. In addition, due to the attractive 
interaction between atoms, the density of the central peak would increase rapidly. In this situation, it is possible that 
the condensates would collapse and even explode in a subsequent time, in analogous with the dynamic process of the 
collapsing and exploding atoms jl5| by switching the interaction from repulsive to attractive. 
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FIGURE CAPTION 



Fig. 1 Displayed is the density of the condensate at x = vs time t, after the optical lattice is switched off. Here 
the density n(x — 0,t) is in units of NA^. We can see clearly in the figure that there is a decay and revival of the 
density oscillation. In addition, there is a periodicity of the density due to the confinement of the magnetic trap. 

Fig. 2 (a)-(d) show the evolution of the density distribution with time t, after the optical lattices are switched off. 
The density distributions are shown at t = 0, 0.1%/u) x , 0.3tt/ui x , and 0.5tt/uj x . The emergence and motion of the 
n = ±1 peaks are clearly shown in these figures. Here the density distribution n(x, t) is in units of NA^, while the 
location x is in units of d, i.e. the distance between two neighboring condensates. 

Fig. 3 Displayed is the motion of the n = 1 peak, after the optical lattices are switched off. Here the location 
x (t) n=1 of n = 1 peak is in units of d. The solid line is the result calculated from the classical harmonic motion given 
by Eq. (|TT|). The squares show the motion of n = 1 peak obtained from the numerical result given by Eq. (§). We 
see that the classical harmonic motion agrees quite well with the numerical result. 

Fig. 4 Displayed is the density of the condensate at x = kudfl vs time t, after the optical lattices are switched off. 
Here the density n(x, t) is in units of NA^. There is a phenomenon of the decay and revival of the density oscillation. 

Fig. 5 Displayed is the time of the disappearance of the density oscillation for different locations. Here the location 
x is in units of d, while the time t is in units of the time scale t w . The solid line is obtained from the formula (pC)), 
while the squares show the result obtained directly from the numerical result given by Eq. (^). We see that the 
analytical formula ( pfj| ) can give a well description for the disappearance of the density oscillation. 

Fig. 6 Displayed is the ratio between n (x = 0,t m , s) and n (x — 0, t m , s = 5) for the interacting and non- interacting 
models. 
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